Genomic trajectories of a near-extinction event in the Chatham Island black robin.

BACKGROUND: Understanding the micro--evolutionary response of populations to demographic declines is a major goal in evolutionary and conservation biology. In small populations, genetic drift can lead to an accumulation of deleterious mutations, which will increase the risk of extinction. However, demographic recovery can still occur after extreme declines, suggesting that natural selection may purge deleterious mutations, even in extremely small populations. The Chatham Island black robin (Petroica traversi) is arguably the most inbred bird species in the world. It avoided imminent extinction in the early 1980s and after a remarkable recovery from a single pair, a second population was established and the two extant populations have evolved in complete isolation since then. Here, we analysed 52 modern and historical genomes to examine the genomic consequences of this extreme bottleneck and the subsequent translocation. RESULTS: We found evidence for two-fold decline in heterozygosity and three- to four-fold increase in inbreeding in modern genomes. Moreover, there was partial support for temporal reduction in total load for detrimental variation. In contrast, compared to historical genomes, modern genomes showed a significantly higher realised load, reflecting the temporal increase in inbreeding. Furthermore, the translocation induced only small changes in the frequency of deleterious alleles, with the majority of detrimental variation being shared between the two populations. CONCLUSION: Our results highlight the dynamics of mutational load in a species that recovered from the brink of extinction, and show rather limited temporal changes in mutational load. We hypothesise that ancestral purging may have been facilitated by population fragmentation and isolation on several islands for thousands of generations and may have already reduced much of the highly deleterious load well before human arrival and introduction of pests to the archipelago. The majority of fixed deleterious variation was shared between the modern populations, but translocation of individuals with low mutational load could possibly mitigate further fixation of high-frequency deleterious variation.

Evolution of gynaecologists' practices regarding the implementation of Swiss legislation on maternity protection at work between 2008 and 2017.

BACKGROUND: In accordance with the International Labour Organization’s Maternity Protection Convention (No. 183) and European Union Directive 92/857CEE (1992), Switzerland’s Labour Law and its Maternity Protection Ordinance (OProMa) aim to protect the health of pregnant employees and their future children while enabling them to pursue their working activities. Gynaecologists-obstetricians have a key role in this legislation, particularly through the prescription of preventive leave for patients who would otherwise face dangerous or arduous tasks in the absence of an adequate risk analysis or suitable protective measures. However, international and national literature suggests that gynaecologists-obstetricians may encounter difficulties in fulfilling their role. AIMS: This study aimed to: (1) describe the practices and difficulties encountered by gynaecologists-obstetricians in the practical implementation of the OProMa; and (2) compare the evolution of these practices and difficulties between 2008 and 2017. METHODS: A survey by questionnaire was conducted in 2008 and repeated in 2017. Both surveys focused on gynaecologists-obstetricians working in the French-speaking part of Switzerland (in private practices, hospitals or both). Descriptive and comparative analyses were carried out. RESULTS: 83 gynaecologists-obstetricians responded in 2008 and 93 in 2017: response rates of 47% and 32%, respectively. In 2017, gynaecologists-obstetricians were more likely to ask questions about occupational risks faced by their patients when consulted by working mothers about their pregnancies. The estimated percentage of patients exposed to an occupational risk remained constant (20% in 2008 and 22% in 2017). Communication and collaboration with employers were reported to be difficult in both surveys, even though these are key elements in the implementation of the OProMa. Collaboration with occupational physicians, however, was more frequent in 2017. CONCLUSION: In 2017, gynaecologists-obstetricians showed a greater awareness of occupational risks and collaborated more frequently with occupational health specialists. However, the application of the OProMa remained limited over the studied time period. Improving training of gynaecologists-obstetricians in this field could be a significant factor in encouraging better implementation of the current legislation. Moreover, gynaecologists-obstetricians need to be given the necessary support to enable their clinical practice to evolve towards a more preventive type of medicine. Collaboration with relevant stakeholders, including occupational physicians, midwives and workers, should be encouraged.

Conservation genetics of two threatened frogs from the Mambilla highlands, Nigeria.

Amphibians are the vertebrate group with the highest number of species threatened with extinction, and habitat loss and fragmentation are considered to be among the leading causes of their declines and extinctions. Little is known of the population biology of amphibian species inhabiting montane forests in Central and West Africa, where anthropogenic activities such as farming and cattle raising are major threats to native biodiversity. We used Amplified Fragment Length Polymorphisms (AFLPs) to assess the population genetic structure of two poorly known species, Cardioglossa schioetzi and Leptodactylodon bicolor (both in the Arthroleptidae), in and around Ngel Nyaki Forest Reserve on the Mambilla Plateau in eastern Nigeria. The landscape comprises continuous forest on steep slopes and small riparian forest fragments in a grassland matrix. While increased fragmentation is well documented for these and other forests in the mountains of Cameroon and Nigeria over the past century, there are no previous assessments of the impact of forest fragmentation on montane amphibian populations in this region. Our estimates of genetic diversity are similar across populations within each species with levels of heterozygosity values consistent with local population declines. Except for a pair of populations (C. schioetzi) we did not observe genetic differentiation between forest and riparian forest fragment populations, nor across sites within continuous forest (L. bicolor). Our results demonstrate recent gene flow between forest fragments and the adjacent protected forests and suggest that small forest corridors connecting these may lessen the genetic consequences of at least 30 years of intense and severe fragmentation in Ngel Nyaki.

Isolation of polymorphic microsatellite loci in the New Zealand endemic sand-binder, Ficinia spiralis (Cyperaceae).

PREMISE OF THE STUDY: Ficinia spiralis (Cyperaceae) is a declining sand-binding sedge of ecological and cultural importance. Microsatellite primers were developed in F. spiralis to investigate how population genetic structure is related to the pronounced morphological, physiological, and ecological variation observed in this species. METHODS AND RESULTS: A 454 shotgun-sequencing approach was used to generate 157,274 raw sequence reads, 536 of which contained microsatellites. After initial primer testing for 40 loci, 14 polymorphic loci were isolated, containing five to 27 alleles per locus. Ten of these loci also amplified in a congener, F. nodosa. CONCLUSIONS: These loci will enable the assessment of the population genetic structure of F. spiralis, improving our understanding of the population processes underlying the observed morphological, physiological, and ecological variation in this endemic species. As the first microsatellites developed in Ficinia, these loci are a valuable resource for population genetic studies within this genus.

A comparative investigation of azole susceptibility in Candida isolates from vulvovaginal candidiasis and recurrent vulvovaginal candidiasis patients in Ghana.

Vulvovaginal candidiasis (VVC) and recurrent vulvovaginal candidiasis (RVVC) affect millions of women and are typically treated with azoles. We know little about azole susceptibility of Candida species from VVC versus RVVC patients, and nothing about African isolates. We have investigated the susceptibility of Candida isolates from Ghana to fluconazole, itraconazole and/or voriconazole. The percentage of Candida albicans isolates showing susceptibility was significantly lower in RVVC than VVC patients. Isolates of Candida parapsilosis and Candida tropicalis showed a similar trend. For Candida glabrata there was no observed difference. The data indicate a decreased susceptibility in selected Candida species from RVVC patients.

Explaining large mitochondrial sequence differences within a population sample.

Mitochondrial DNA sequence is frequently used to infer species' boundaries, as divergence is relatively rapid when populations are reproductively isolated. However, the shared history of a non-recombining gene naturally leads to correlation of pairwise differences, resulting in mtDNA clusters that might be mistaken for evidence of multiple species. There are four distinct processes that can explain high levels of mtDNA sequence difference within a single sample. Here, we examine one case in detail as an exemplar to distinguish among competing hypotheses. Within our sample of tree weta (Hemideina crassidens; Orthoptera), we found multiple mtDNA haplotypes for a protein-coding region (cytb/ND1) that differed by a maximum of 7.9%. From sequencing the whole mitochondrial genome of two representative individuals, we found evidence of constraining selection. Heterozygotes were as common as expected under random mating at five nuclear loci. Morphological traits and nuclear markers did not resolve the mtDNA groupings of individuals. We concluded that the large differences found among our sample of mtDNA sequences were simply owing to a large population size over an extended period of time allowing an equilibrium between mutation and drift to retain a great deal of genetic diversity within a single species.

Building strong relationships between conservation genetics and primary industry leads to mutually beneficial genomic advances.

Several reviews in the past decade have heralded the benefits of embracing high-throughput sequencing technologies to inform conservation policy and the management of threatened species, but few have offered practical advice on how to expedite the transition from conservation genetics to conservation genomics. Here, we argue that an effective and efficient way to navigate this transition is to capitalize on emerging synergies between conservation genetics and primary industry (e.g., agriculture, fisheries, forestry and horticulture). Here, we demonstrate how building strong relationships between conservation geneticists and primary industry scientists is leading to mutually-beneficial outcomes for both disciplines. Based on our collective experience as collaborative New Zealand-based scientists, we also provide insight for forging these cross-sector relationships.

The isolation of Candida rugosa and Candida mesorugosa from clinical samples in Ghana.

Members of the Candida rugosa species complex have been described as emerging fungal pathogens and are responsible for a growing number of Candida infections. In this communication we report the isolation of Candida rugosa and Candida mesorugosa in Ghana. To the best of our knowledge this is the first description of this species complex from a clinical setting in Africa.The isolates were identified on the basis of their rRNA gene internal transcribed spacer (ITS) sequences. For one isolate, obtained from sputum, the sequence grouped well with that of C. rugosa. Two other isolates from urine had sequences that grouped with Candida mesorugosa. Morphologically, C. rugosa formed white, wrinkled, and flat colonies on Sabouraud Dextrose Agar (SDA), whereas C. mesorugosa formed white, smooth colonies. On chromogenic medium, the isolates formed small, dry greenish-blue colonies with a pale or white border, similar to C. albicans. The C. rugosa isolate produced pseudohyphae in human serum and on CMA-Tween 80 agar. In contrast, the C. mesorugosa isolates did not generate pseudohyphae in human serum, but generated a few pseudohyphae with abundant blastoconidia on CMA-Tween 80 agar. Growth was observed at 37 °C and 42 °C but not at 45 °C.The two C. mesorugosa isolates had Minimum Inhibitory Concentrations (MICs) of 6 and 48 µg ml(-1) for fluconazole and are thus resistant. The C. rugosa isolate had an MIC of 24 µg ml(-1), indicative of resistance. All three isolates were susceptible to itraconazole and voriconazole (with respective MICs of < 0.125 µg ml(-1)).

Population Structure of Candida albicans from Three Teaching Hospitals in Ghana.

Previous studies on Candida species in a clinical setting in Ghana have shown a prevalence of Candida albicans. Despite this, very little is known about the various strain types and their population genetic structure. In this study three microsatellite loci, CAI, CAIII and CAVI, were used to investigate the population genetic structure of C. albicans from clinical isolates in Ghana. In all, 240 clinically unrelated C. albicans isolates were recovered from patients reporting at three teaching hospitals. All the isolates were heterozygous for at least one of the three loci, except for one isolate, which was homozygous for all three loci. Sixty-seven unique alleles and 240 different genotypes were generated by the three polymorphic microsatellite loci, resulting in a very high discriminatory potential of approximately 0.98. There was no significant difference in allele frequencies from the small number of anatomical sites sampled, regardless of the host conditions although high genotypic diversities were observed among the isolates. There was evidence for clonal reproduction, including over-expression of observed heterozygotes across the populations. The populations deviated significantly from Hardy-Weinberg equilibrium and pair-wise genotypic linkage disequilibria comparisons across the three loci were significant, also suggesting a clonal population. The overall Wright FIS for the three loci was negative, and the overall FST value was not significantly different from zero for the three loci analyzed, indicating a clonal and homogeneous population across the three sampling locations from Ghana.

Extinct New Zealand megafauna were not in decline before human colonization.

The extinction of New Zealand's moa (Aves: Dinornithiformes) followed the arrival of humans in the late 13th century and was the final event of the prehistoric Late Quaternary megafauna extinctions. Determining the state of the moa populations in the pre-extinction period is fundamental to understanding the causes of the event. We sampled 281 moa individuals and combined radiocarbon dating with ancient DNA analyses to help resolve the extinction debate and gain insights into moa biology. The samples, which were predominantly from the last 4,000 years preceding the extinction, represent four sympatric moa species excavated from five adjacent fossil deposits. We characterized the moa assemblage using mitochondrial DNA and nuclear microsatellite markers developed specifically for moa. Although genetic diversity differed significantly among the four species, we found that the millennia preceding the extinction were characterized by a remarkable degree of genetic stability in all species, with no loss of heterozygosity and no shifts in allele frequencies over time. The extinction event itself was too rapid to be manifested in the moa gene pools. Contradicting previous claims of a decline in moa before Polynesian settlement in New Zealand, our findings indicate that the populations were large and stable before suddenly disappearing. This interpretation is supported by approximate Bayesian computation analyses. Our analyses consolidate the disappearance of moa as the most rapid, human-facilitated megafauna extinction documented to date.

Human-assisted spread of a maladaptive behavior in a critically endangered bird.

Conservation management often focuses on counteracting the adverse effects of human activities on threatened populations. However, conservation measures may unintentionally relax selection by allowing the 'survival of the not-so-fit', increasing the risk of fixation of maladaptive traits. Here, we report such a case in the critically-endangered Chatham Island black robin (Petroica traversi) which, in 1980, was reduced to a single breeding pair. Following this bottleneck, some females were observed to lay eggs on the rims of their nests. Rim eggs left in place always failed to hatch. To expedite population recovery, rim eggs were repositioned inside nests, yielding viable hatchlings. Repositioning resulted in rapid growth of the black robin population, but by 1989 over 50% of all females were laying rim eggs. We used an exceptional, species-wide pedigree to consider both recessive and dominant models of inheritance over all plausible founder genotype combinations at a biallelic and possibly sex-linked locus. The pattern of rim laying is best fitted as an autosomal dominant Mendelian trait. Using a phenotype permutation test we could also reject the null hypothesis of non-heritability for this trait in favour of our best-fitting model of heritability. Data collected after intervention ceased shows that the frequency of rim laying has strongly declined, and that this trait is maladaptive. This episode yields an important lesson for conservation biology: fixation of maladaptive traits could render small threatened populations completely dependent on humans for reproduction, irreversibly compromising the long term viability of populations humanity seeks to conserve.

The half-life of DNA in bone: measuring decay kinetics in 158 dated fossils.

Claims of extreme survival of DNA have emphasized the need for reliable models of DNA degradation through time. By analysing mitochondrial DNA (mtDNA) from 158 radiocarbon-dated bones of the extinct New Zealand moa, we confirm empirically a long-hypothesized exponential decay relationship. The average DNA half-life within this geographically constrained fossil assemblage was estimated to be 521 years for a 242 bp mtDNA sequence, corresponding to a per nucleotide fragmentation rate (k) of 5.50 × 10(-6) per year. With an effective burial temperature of 13.1°C, the rate is almost 400 times slower than predicted from published kinetic data of in vitro DNA depurination at pH 5. Although best described by an exponential model (R(2) = 0.39), considerable sample-to-sample variance in DNA preservation could not be accounted for by geologic age. This variation likely derives from differences in taphonomy and bone diagenesis, which have confounded previous, less spatially constrained attempts to study DNA decay kinetics. Lastly, by calculating DNA fragmentation rates on Illumina HiSeq data, we show that nuclear DNA has degraded at least twice as fast as mtDNA. These results provide a baseline for predicting long-term DNA survival in bone.

Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies.

One of the most common questions asked before starting a new population genetic study using microsatellite allele frequencies is "how many individuals do I need to sample from each population?" This question has previously been answered by addressing how many individuals are needed to detect all of the alleles present in a population (i.e. rarefaction based analyses). However, we argue that obtaining accurate allele frequencies and accurate estimates of diversity are much more important than detecting all of the alleles, given that very rare alleles (i.e. new mutations) are not very informative for assessing genetic diversity within a population or genetic structure among populations. Here we present a comparison of allele frequencies, expected heterozygosities and genetic distances between real and simulated populations by randomly subsampling 5-100 individuals from four empirical microsatellite genotype datasets (Formica lugubris, Sciurus vulgaris, Thalassarche melanophris, and Himantopus novaezelandia) to create 100 replicate datasets at each sample size. Despite differences in taxon (two birds, one mammal, one insect), population size, number of loci and polymorphism across loci, the degree of differences between simulated and empirical dataset allele frequencies, expected heterozygosities and pairwise F(ST) values were almost identical among the four datasets at each sample size. Variability in allele frequency and expected heterozygosity among replicates decreased with increasing sample size, but these decreases were minimal above sample sizes of 25 to 30. Therefore, there appears to be little benefit in sampling more than 25 to 30 individuals per population for population genetic studies based on microsatellite allele frequencies.

Profiling the dead: generating microsatellite data from fossil bones of extinct megafauna--protocols, problems, and prospects.

We present the first set of microsatellite markers developed exclusively for an extinct taxon. Microsatellite data have been analysed in thousands of genetic studies on extant species but the technology can be problematic when applied to low copy number (LCN) DNA. It is therefore rarely used on substrates more than a few decades old. Now, with the primers and protocols presented here, microsatellite markers are available to study the extinct New Zealand moa (Aves: Dinornithiformes) and, as with single nucleotide polymorphism (SNP) technology, the markers represent a means by which the field of ancient DNA can (preservation allowing) move on from its reliance on mitochondrial DNA. Candidate markers were identified using high throughput sequencing technology (GS-FLX) on DNA extracted from fossil moa bone and eggshell. From the 'shotgun' reads, >60 primer pairs were designed and tested on DNA from bones of the South Island giant moa (Dinornis robustus). Six polymorphic loci were characterised and used to assess measures of genetic diversity. Because of low template numbers, typical of ancient DNA, allelic dropout was observed in 36-70% of the PCR reactions at each microsatellite marker. However, a comprehensive survey of allelic dropout, combined with supporting quantitative PCR data, allowed us to establish a set of criteria that maximised data fidelity. Finally, we demonstrated the viability of the primers and the protocols, by compiling a full Dinornis microsatellite dataset representing fossils of c. 600-5000 years of age. A multi-locus genotype was obtained from 74 individuals (84% success rate), and the data showed no signs of being compromised by allelic dropout. The methodology presented here provides a framework by which to generate and evaluate microsatellite data from samples of much greater antiquity than attempted before, and opens new opportunities for ancient DNA research.

Genetic analyses reveal hybridization but no hybrid swarm in one of the world's rarest birds.

Hybridization facilitated by human activities has dramatically altered the evolutionary trajectories of threatened taxa around the globe. Whereas introduced mammalian predators and widespread habitat loss and degradation clearly imperil the recovery and survival of the New Zealand endemic black stilt or kaki (Himantopus novaezelandiae), the risk associated with hybridization between this critically endangered endemic and its self-introduced congener, the pied stilt or poaka (Himantopus himantopus leucocephalus) is less clear. Here, we combine Bayesian admixture analyses of microsatellite data with mitochondrial DNA sequence data to assess the levels of hybridization and introgression between kaki and poaka. We show that birds classified as hybrids on the basis of adult plumage are indeed of hybrid origin and that hybridization between kaki and poaka is both extensive and bidirectional. Despite this, we found almost no evidence for introgression from poaka to kaki, thus negating the popular belief that kaki represent a hybrid swarm. To our knowledge, ours represents the first comprehensive study to document a lack of widespread introgression for a species at risk despite a recent history of extensive bidirectional human-induced hybridization. We attribute this rather surprising result, in part, to reduced reproductive success in female hybrids combined with a transient male-biased kaki sex ratio. To maximize the evolutionary potential of kaki, we use these data to recommend conservation management activities aimed to maintain the genetic integrity and to maximize the genetic diversity of this iconic rare bird.

Merging ancient and modern DNA: extinct seabird taxon rediscovered in the North Tasman Sea.

Ancient DNA has revolutionized the way in which evolutionary biologists research both extinct and extant taxa, from the inference of evolutionary history to the resolution of taxonomy. Here, we present, to our knowledge, the first study to report the rediscovery of an 'extinct' avian taxon, the Tasman booby (Sula tasmani), using classical palaeontological data combined with ancient and modern DNA data. Contrary to earlier work, we show an overlap in size between fossil and modern birds in the North Tasman Sea (classified currently as S. tasmani and Sula dactylatra fullagari, respectively). In addition, we show that Holocene fossil birds have mitochondrial control region sequences that are identical to those found in modern birds. These results indicate that the Tasman booby is not an extinct taxon: S. dactylatra fullagari O'Brien & Davies, 1990 is therefore a junior synonym of Sula tasmani van Tets, Meredith, Fullagar & Davidson, 1988 and all North Tasman Sea boobies should be known as S. d. tasmani. In addition to reporting the rediscovery of an extinct avian taxon, our study highlights the need for researchers to be cognizant of multidisciplinary approaches to understanding taxonomy and past biodiversity.

Identification of microsatellites from an extinct moa species using high-throughput (454) sequence data.

Genetic variation in microsatellites is rarely examined in the field of ancient DNA (aDNA) due to the low quantity of nuclear DNA in the fossil record together with the lack of characterized nuclear markers in extinct species. 454 sequencing platforms provide a new high-throughput technology capable of generating up to 1 gigabases per run as short (200-400-bp) read lengths. 454 data were generated from the fossil bone of an extinct New Zealand moa (Aves: Dinornithiformes). We identified numerous short tandem repeat (STR) motifs, and here present the successful isolation and characterization of one polymorphic microsatellite (Moa_MS2). Primers designed to flank this locus amplified all three moa species tested here. The presented method proved to be a fast and efficient way of identifying microsatellite markers in ancient DNA templates and, depending on biomolecule preservation, has the potential of enabling high-resolution population genetic studies of extinct taxa. As sequence read lengths of the 454 platforms and its competitors (e.g., the SOLEXA and SOLiD platforms) increase, this approach will become increasingly powerful in identifying microsatellites in extinct (and extant) organisms, and will afford new opportunities to study past biodiversity and extinction processes.

Characterization of microsatellite markers for Sycoscapter nonpollinating fig wasps.

We isolated 18 microsatellites from Sycoscapter australis, a nonpollinating fig wasp that develops in figs of Ficus macrophylla, and assessed their variability in 20 wasps. We further optimized nine of these loci for use in three other Sycoscapter species that develop in Ficus rubiginosa figs and assessed their variability in 47-140 wasps per species. These are the first microsatellites developed for nonpollinating fig wasps and show sufficient polymorphism to become important tools in evolutionary and genetical studies of Sycoscapter wasps.

Development of polymorphic microsatellite markers for the New Zealand black stilt (Himantopus novaezelandiae) and cross-amplification in the pied stilt (Himantopus himantopus leucocephalus).

Eight polymorphic microsatellite primer pairs were developed for the critically endangered New Zealand black stilt, Himantopus novaezelandiae, representing the first microsatellite markers available for birds in the family Recurvirostridae. The number of alleles ranged from two to four per locus. Observed and expected heterozygosities ranged from 0.30 to 0.80 and from 0.37 to 0.70, respectively. All eight loci were polymorphic in the related species Himantopus himantopus leucocephalus, indicating these primer pairs may be useful for additional taxa in the globally distributed genus Himantopus.

Causes of size homoplasy among chloroplast microsatellites in closely related Clusia species.

Chloroplast DNA sequences and microsatellites are useful tools for phylogenetic as well as population genetic analyses of plants. Chloroplast microsatellites tend to be less variable than nuclear microsatellites and therefore they may not be as powerful as nuclear microsatellites for within-species population analysis. However, chloroplast microsatellites may be useful for phylogenetic analysis between closely related taxa when more conventional loci, such as ITS or chloroplast sequence data, are not variable enough to resolve phylogenetic relationships in all clades. To determine the limits of chloroplast microsatellites as tools in phylogenetic analyses, we need to understand their evolution. Thus, we examined and compared phylogenetic relationships of species within the genus Clusia, using both chloroplast sequence data and variation at seven chloroplast microsatellite loci. Neither ITS nor chloroplast sequences were variable enough to resolve relationships within some sections of the genus, yet chloroplast microsatellite loci were too variable to provide any useful phylogenetic information. Size homoplasy was apparent, caused by base substitutions within the microsatellite, base substitutions in the flanking regions, indels in the flanking regions, multiple microsatellites within a fragment, and forward/reverse mutations of repeat length resulting in microsatellites of identical base composition that were not identical by descent.